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ABSTRACT 


A technique has been developed which makes it possible to 
study (n,n'y) reactions using very small scattering samples in 
a close geometry. To show that reliable y-ray yield curves can 
be obtained with this method, the excitation curve for the 847 
keV y-ray from °°Fe has been measured. The resulting yield curve 
is in good agreement with previous measurements obtained by 
conventional techniques. 

The method has been applied to the study of levels above 
900 keV excitation in !5*Sm and !©°Gd populated by the (n,n'y) 
reaction. De-excitation y-rays were observed at several neutron 
energies between 0.9 and 2.3 MeV. The y-ray production cross 
sections have been obtained at these neutron energies and were 
used to determine neutron inelastic scattering cross sections for 
reactions leading to several states in these nuclei. The experi- 
mental results were compared with the predictions of the statis- 
tical theory to establish level spins. 

The 1°*Sm nucleus: Previously known spin assignments have 
been confirmed: 1 for the 921 keV level, 3 for the 1012 keV 
level, 0° for the 1099 keV level, 1” for the 1475 keV level and 
3 for the 1585 keV level. The 1178 and 1182 keV levels were 
clearly resolved and assigned spins of ai and 5 , respectively. 


Existence of levels at 1202, 1286. 1338, 1515, 1539, 1707 and 
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1755 keV has been confirmed. A level at 1674 keV is tentatively 
proposed. Spin assignments of 4* and au have been made for the 
1338 and 1539 keV levels, respectively. Tentative assignments 
are: 0° to the 1202 keV level, 37 to the 1286 keV level, 

3 to the 1515 keV level, O* to the 1674 keV level, 4° to 
the 1707 keV level and 2” to the 1755 keV level. 

The 1®°Gd nucleus: Spin assignments of 2° and 3° were 
obtained for the 989 and 1290 keV levels, respectively, in 
agreement with previous results. An assignment of 3’ and 10 has 
been made for the 1058 and 1225 keV states, respectively. The 
levels at 1071 and 1149 keV were tentatively assigned a J" = 4". 
Seven new levels are proposed at 1262, 1351, 1378, 1381, 1569, 
1587 and 1599 keV. A spin of 1. was assigned to the 1351 keV 
level. Tentative assignments are: 1 to the 1262 keV level, 0 
or 1” to the 1378 and 1381 keV levels, and 1” or 2° to the 1569, 
1587 and 1599 keV levels. 


Reduced transition probability ratios for y-rays depopulating 


members of rotational bands observed in !°*Sm and 1°°Gd were 


compared to the predictions of the Bohr-Mottelson collective model. 


The results indicate that these nuclei are good rotors. 
The ’Li(p,n)’Be reaction was used as a neutron source. 


Analysis of the y-ray data requires an accurate knowledge of the 


angular distribution of neutrons from this reaction. Consequently 


a systematic measurement of absolute differential cross section as 


a function of both energy and angle for the ’Li(p,n,)’Be and 
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’Li(p.ni)’Be reactions has been carried out at proton energies from 


2.60 to 5.40 MeV and laboratory angles between 0° and 150°. 
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CHAPTER I 


INTRODUCTION 


Measurements of y-ray energies and intensities emitted after 
neutron inelastic scattering, coupled with the application of the 
Statistical theory of nuclear reactions, can provide useful infor- 
mation about the energy, decay mode and spin of excited states of 
the scattering nuclei. 

The use of neutrons to excite levels for decay studies is 
experimentally favorable over that of charged particles. With 
neutrons one can, by adjusting their energy, excite and observe 
each level in turn so that the y-spectrum is not complicated by 
radiation due to higher levels. With charged particles, however, 
higher bombarding energies are required to overcome the Coulomb 
barrier, resulting in excitation of a large number of levels at 
Once and increasing the probability of occurance of competing 
reactions which contribute to the complexities of the y-spectrum. 

The possibility of exciting only a few final states in 
neutron inelastic scattering leads to a simple situation in which 
the theory can predict the magnitude of the cross section within 
10% (Vo 68). This remarkable success of the theory has greatly 
enhanced the effectiveness of these reactions for establishing level 
Spins. 


The number of nuclei which can be studied by neutron inelastic 
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scattering is normally limited by the necessity of using large 
scattering samples containing 1 -2 mol of metayi atl Many of the 
elements of interest, especially in the A > 100 mass region where 
the statistical theory should be particularly applicable, occur 
in an isotopically impure form and this obviously complicates the 
interpretation of results. It is often prohibitively expensive 
or physically impossible to obtain such large quantities of the 
enriched isotope. 

This thesis describes a new method which has been developed 
to carry out (n,n'y) reactions with only 0.01 mol of scattering 
material. This has been accomplished by placing the scatterer in 
close proximity to the neutron source to make the neutron flux as 
large as possible in the scatterer. Two difficulties then occur. 
Firstly, the neutron beam through the scatterer is neither mono- 
energetic nor of constant flux; this is because energy and cross 
section are quite angle dependent for all the commonly used neutron 
producing reactions (Go 63). Secondly, there is no possibility of 
shielding the y-detector from the neutron beam. An analysis routine 
has been developed which overcomes the first difficulty. In 
addition, time-of-flight gating of the y-ray spectra has proven to 
be adequate to reduce background to acceptable levels despite the 
lack of shielding. In order to check whether this method would give 
reliable cross section information, the excitation curve of the 
847 keV y-ray of °°Fe was measured using this technique with a 


3.147 g sample of natural iron. When the results were compared 
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with those of the conventional technique, good agreement was 
obtained. 

The method was then applied to the study of nuclei in the 
rare-earth region, namely !°*Sm and 1°°Gd. Although the structure 
of '°*Sm has been the subject - several studies (Ke 64, Yo 65, 

Bj 66, Ze 66, Se 66; Ve 68, Da 71), relatively little information 
about the levels of !®°Gd is available (Yo 65, Sh 66, Bl 67). 

Prior to the initiation of the present work, knowledge of the 

decay properties of both nuclei were incomplete and severely limited 
by the poor resOlution of Nal detectors. Moreover, level spins were 
not well established. The reactions previously used to study these 
nuclei preferentially excite some states depending on their structure. 
The present investigation has been carried out using the (n,n'y) 
reaction which has the advantage of populating levels independent 

of their character, provided that the statistical theory require- 
ments are satisfied. Furthermore, the measured inelastic cross 
sections can be compared with the predictions of the statistical 
theory to determine level spins. 

In the present study we have utilized the high resolution of 
Ge(Li) detectors to observe y-rays following neutron-inelastic 
scattering on !°*Sm and '°°Gd. Level excitation energies and decay 
modes have been inferred from measurements of y-ray energies and 
intensities, respectively. Inelastic cross sections for the observed 
levels were determined at incident neutron energies of 1 to 2 MeV 


and compared with theoretical calculations to obtain level spins. 
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Reduced y-ray transition probabilities have been interpreted in the 
framework of the Bohr and Mottelson (Bo 53) collective model of 
deformed nuclei. 

It is worth mentioning that while the analysis of the present 
data was in progress, the y-decay work of D'Auria et al. (Da 71) 
came to our attention prior to its publication and the results were 
exchanged. 

An outline of the statistical theory of nuclear reactions is 
given in Chapter II. The small sample method has been described in 
Chapter III, where a measurement of the excitation curve of the 
847 keV y-ray from °°Fe and its comparison to previous results are 
reported. Chapter IV contains the application of this method to 
the study of levels in ?°*Sm and '®°Gd, whereas the conclusions of 


the present work are summarized in Chapter V. 
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THE STATISTICAL THEORY OF NUCLEAR REACTIONS 


2.1 Introduction 


The statistical theory is based on the compound nucleus (CN) 
picture of nuclear reactions. According to this picture the incident 
projectile is absorbed by the target nucleus to form an excited 
compound system. This compound nucleus subsequently decays, generally 
after a time long compared with the nuclear transit time, into the 
energetically allowed channels. At low incident energies, compound 
nucleus processes are characterized by the appearance of isolated 
resonances in their excitation functions; these resonances correspond 
to levels in the CN. As the energy increases, the number of reso- 
nances increases and their widths [ become greater than the spacings 
D between them so that they overlap with each other. If the energy 
Spread of the incident beam AE is such that AE >> I, then the 
measured cross sections represent an energy average over many 
resonances that are covered by the energy spread of the incident 
beam. The statistical theory of nuclear reactions deals with such 
average cross sections. 

aaenemses aa development of the statistical theory is given 
in numerous publications which are referred to in the recent review 
article by Vogt (Vo 68). In this chapter, we give a brief outline 


of the ideas involved in the statistical model treatment. 
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2.2 Statistical Theory Predictions for Average Cross Sections 
Consider the reaction 
d.> 6 


where a labels the incident particle, the target nucleus and 
their state of excitation. a is called the entrance 
channel. 

8 labels the outgoing particle, the residual nucleus and 
their state of excitation. 8 is one of many possible 
exit channels. 

In the absence of spins, the cross section for this reaction can 


be written as (Fr 55) 
tggle Te 
Oo = | Sug 7 Sug | 20] 


where Ke is the wave number of the relative motion of the pair 


Q, and Sy is the scattering matrix elemen= which is the 


B 
amplitude of the outgoing wave in the channel 8 when there is 

an incoming wave of unit amplitude in the channel a. In 

general, this reaction can occur directly without passing through 
the CN state and/or by forming a CN which decays into the channel 
8. Thus the cross section Fug will have two components; a direct 
interaction (DI) component and a CN component. The DI component 
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varies smoothly with energy. If the CN is formed at a high excitation 
where the levels strongly overlap then the CN component of O48 will 
fluctuate rapidly with energy. These fluctuations are indicative 

of rapid varictions of S08" It is then convenient to decompose the 


S-matrix into two parts -- S which varies smoothly with energy and 


S which fluctuates rapidly with energy. Thus 
S es. +S a2] 
Eq. (2.1) then becomes 


ive eae Se : 
oe ae | Se Sug Sia! i2ean 
Ol 
If the energy spread of the incident beam is large enough to 
excite many levels of the CN, then the measured cross section 
represents an energy average. Therefore we average expression (2.3) 


over the beam energy spread to get: 


ive emule ne is 

<0 5a = kee <I (Sug Sg) ae > [2.4] 
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This hypothesis allows the factorization of the cross section into 


a formation part and a decay part: 


CN 


<o 
ap 


Sr) ie Oy Pa 


where Oy is the cross section for formation of the CN and P, is 


B 
its decay probability through the channel 8. 


Since the CN must decay, then 


and hence 


But from (2.8) and (2.9) we have 


CN. es tees 
<O ya > wr ke <|Syel? 
therefore 
PRs i P 
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From conservation of flux, it follows that 


the S-martix ts unitary, 1.e. 


i 
Sap Sug 1 


[210] 
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The corresponding energy-averaged expression is 


y <Sia eae (ed 
8 
then we get 
<tr erer ee : 
1 Sg * Sop) (S.g* Sy)? = ] 
J [Sygl? #2 <ISygl* = 1 [2.16] 
B 
Substituting (2.16) in (2.14) we obtain 
x 7 
o, 5 Ke ft LDagh ) [7] 
a 


This equation shows that the elements of the DI part of the S- 
matrix are capable of determining the CN formation cross section. 
In pursuing the consequences of this striking observation, we 
note that eq. (2.17) can be written in terms of a transmission 


coefficient Te : 


| 
rT 


ie 1 Sag [2.18] 


e) 


which expresses the loss so far as the DI processes are 


concerned. Eq. (2.17) then takes the form 


o = pT [2.19] 
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The optical model permits, in a phenomenological way, 
the calculation of the diagonal elements Sy of the S-matrix. 
These are related to the optical model phase shifts 84 by the 


relation 


The transmission coefficients (2.18) are then approximated 


by the optical model transmission coefficient 
; 2. 
T = 1- aa 


For low energy reactions where nonelastic direct reactions 
(a # 8) are negligible, eq. (2.21) is considered to be a good 
approximation. 


The cross section <oug is now given by: 


To calculate Pe we assume that 


the S-matrix is invariant under time-reversal, 1.e. 


S08 “ oR 


where reversed arrows denote time-reversed channels. This 


implies that the amplitude for the reaction a > 8 is the same 


[2.20] 
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as the amplitude for the time-reversed reaction 8 > a; this is 
known as the reciprocity theorem. (Although spin projections 
will be reversed, we are going to ignore that here since it does 
not effect the final result because we always sum over spin 
projections.) The reciprocity property for the fluctuating part 
of the S-matrix leads to the following relation between cross 
sections: 


2 CNS ee 2 CN 
ky “Chokes k. <oge > [2.24] 


which, together with eq. (2.22), gives 


je ce te h2.201 


le Pot a oe. i226. 
from (2.25) and (2.26) we get 


Pe iP T2/Ty 


which allows us to write: 


Pig U4 idol ea 


where a is constant. If we sum eq. (2.27) over all outgoing 


channels 
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we get 
] 
a = —_—_—— 
iT, 
B 


substituting (2.27) in (2.26) yields 


The CN cross section can now be written in terms of 


transmission coefficients as 


T 
CN a ii 
Soe > aa [ew 


This equation gives the energy average CN cross section in the 


absence of spin. 


2.3 Refinements - The Hauser-Feshbach Formula 


Although it is assumed that the processes of formation and 
decay of the CN are independent, the energy, total angular 
momentum and parity must all be conserved. 
tions are taken into account a refined version of eq. (2.30) is 
obtained; it is often referred to as the Hauser-Feshbach formula 
after the authors who first added conservation of total angular 


momentum and parity to the calculation of neutron inelastic 


cross sections (Ha 52). 


When these restric- 
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To incorporate angular momentum and parity in eq. (2.30) we 
shall use a to denote the pair of particles and their state of 


excitation as well as the following set of quantum numbers: 


S = intrinsic spin of the incident particle 

I = intrinsic spin of the target nucleus 

S = channel spin formed by coupling s with I 

& = orbital angular momentum 

J = total angular momentum formed by coupling & with Ss 
My = z-component of J 

II 3 total parity which is the product of tye with the 


intrinsic parities of the pair a. 
In what follows, we shall assume that all reactions are independent 
of M) and that My can be averaged over and ignored. 
Following the same steps given in §2.2 and applying spin and 
parity conservation restrictions, we can write oui explicitly 


as (Vo 68) 
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aa!” * KF dt Tessar) i Suara: aie 


where the unprimed quantities refer to incoming channels of the 
reaction, the primed quantities refer to the outgoing channels, 
and the double primed quantities in the denominator are summed 

over all channels to which the CN can decay. The sums in eq. 


(2.31) run over the following values: Il is positive or negative; 
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J has the values 0, 1, 2, . . . if the CN has an even mass number 
and the values 1/2, 3/2, 5/2, . . . otherwise; S takes on all 
integer values between |J- S| and J + S if the pair a has the same 
parity as Il and over all odd values otherwise. 

It is clear from eq. (2.31), known as the Hauser-Feshbach 
formula, that the transmission coefficients are independent of 
S and J. This is true only if the optical model potential has no 
spin dependent terms. If a spin-orbit coupling term is present 
in the potential, then the transmission coefficients will depend 
on the angular momentum j of the projectile (j = 2+ s). In 
this case it is more convenient to change to a coupling scheme 
where the intrinsic spin of the incident particle s couples to 
its orbital angular momentum & to give a total angular momentum 
j,; which then couples to the intrinsic spin of the nucleus to give 
the total J. This presents no problem since the transmission 


coefficients in the two coupling schemes are related 


T (a) = 2 (2j+1) (28+1) W? (2sdI3js) T, -(a) [2.32] 
. J 
S2J J 
where T (cx) are now the transmission coefficients that depend 
SQJ 


on s and J to be used in place of T, (a) Ofred. 9 (2-34 Weiss the 
Racah coefficient, and 1); are related to the optical model phase 


shifts 69; (a) by 


Theta) ee elem ert): [2.33] 
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Therefore, with the transformation (2.32) we can still use 
eq. (2.31) to get the cross section even if the optical potential 
contains a spin-orbit term. 

The cross section given by eq. (2.31) is, of course, the 
total cross section, i.e., angle integrated. The differential 
cross section may be derived from somewhat more complicated 
calculations involving Racah algebra; the explicit form is given 
in many books and articles (e.g. Vo 68). 

Inspection of eq. (2.31) indicates that if the angular 
momenta and parities involved in the different reaction channels 
are known, then the numerical calculation of the CN cross section 
requires only a knowledge of the respective transmission coeffic- 
ients, and these, as we have seen, are determined by the diagonal 
elements of the direct part of the S-matrix. In practice, the 
parameters of an appropriate optical model are determined from an 
optimum fit to DI elastic scattering data. Once the potential is 
determined, one can solve Schrddinger's equation with such a 


potential to get the phase shifts. 
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2.4 The Width-Fluctuation Correction 


In cases where the CN levels are non-overlapping and only a 
few channels are open, the HF formula is not sufficiently accurate 
mainly because the HF _ theory inherently presupposes constancy 
in the partial widths of the CN levels. This may be brought out 
more clearly by derivation of the formula of the CN cross section 
by averaging over all the CN resonances, using the Briet-Wigner 
formula (La 57). 

If the Briet-Wigner single level formula is integrated over 
the energy region around the resonance, and averaged over the 
resonances of the same spin in the region of excitation covered 
by the energy spread of the incident beam, one gets (neglecting 


interference between neighboring CN levels) (Mo 64) 
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where Dai is the mean level spacing of resonances with angular 
momentum J and parity Il, the total width a in the denominator 
is the sum over all competing partial widths Fie by" Pig? 


and < denotes an average over resonances u with angular 


“u (dtl) 
momentum J and parity II. 

To put this expression into a form comparable to the HF 
formula (eq. 2.31) we use the following expression for the 


transmission coefficient (Mo 64) 
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The HF formula follows from eqs. (2.34), (2.35) if one sets 


ih Tic! SH tea Sieg 
Gia =y guar STio'* [2.36] 
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It is then apparent that the HF formula ignores the effect of 
fluctuations of CN level widths about their average values. This 
cannot be a good assumption even if widths in different channels 
are completely independent, because fluctuations in the numerator 
are correlated to fluctuations in the denominator. Each term in 
the HF formula should be multiplied by a correction factor of the 


form 
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This width fluctuation correction can be evaluated as a function 
of the average values elie? from knowledge of the form of the 
Statistical distribution of the widths about their average values. 
A useful form is the Porter-Thomas distribution (Po 56) 


P(y) dy = (2ny) erie dy [2.38] 


where" 4g ehul adalah For this distribution Woo! can be expressed 


as a definite integral. Writing Hae in terms of T's, instead of 
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T's, one obtains (Mo 61) 
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x Totes er) ; dy [2.39] 


re 


This integral can be evaluated analytically only for a" < 2 but 
if a >2 it can be evaluated numerically. With very few open 
channels Woo! may be as small as one-half fora #a' (i.e. inelastic 
scattering) but as large as three for a = a' (elastic scattering), 
but it rapidly approaches unity when the number of channels becomes 


large (Mo 64a, Vo 68). 
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CHAPTER III 


THE SMALL SAMPLE METHOD FOR INVESTIGATION OF THE (n,n'y) REACTION 


3.1 Experimental Method 


The experimental arrangement is shown in fig. 3.1. Neutrons 
were produced using the ’Li(p,n)’Be reaction. A metallic ’Li 
target on a 0.254 cm thick Ta backing was bombarded with 0.5 ns 
proton pulses from the 6 MV van de Graaff accelerator in conjunction 
with the Mobley magnet bunching system. The average beam intensity 
was 1 uA during the course of the experiment; the beam spot size 
was approximately 0.2 cm in diameter. To reduce background from 
Surrounding materials, the end of the proton beam tube was made of 
nylon with a gold inner lining. In order to ensure that the proton 
beam was always positioned on the same target spot during the experi- 
ment, the beam was passed through a pre-target slit system. It 
consisted of two gold plates independently mounted on two micrometers 
and connected to two microammeters. A negative bias of 300 V was 
applied to each plate to prevent secondary electrons from leaving 
the target. The Mobley magnet current was adjusted during the 
experiment to ensure a balanced minimum current on the plates. 

The scatterer, placed immediately behind the Ta backing, 
consisted of 3.147 g of natural Fe in the form of a disc, 1.584 cm 
in diameter and 0.204 cm thick. 


Gamma rays were detected in a 48 cm® Ge(Li) detector placed 
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Figure 3.1] Schematic diagram of the experimental set-up. 
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40 cm from the scatterer at an angle of 90° with respect to the 
incident proton beam. Since the detector itself is sensitive to 
fast neutrons, re necessary to discriminate against background 
Y-rays resulting from neutron interactions with the Ge(Li) crystal 
and the surrounding materials. The technique adopted utilizes the 
timing capability of the detector. 

Signals from the Ge(Li) detector were used to obtain both 
energy and timing information. Timing signals were extracted using 
a timing filter amplifier (TFA) and a constant-fraction timing dis- 
criminator (CFTD) and provided the start signal for the time-to- 
amplitude converter (TAC). The stop signal was derived from a 
cylindrical capacitive time pick-off placed in the beam tube 
immediately ahead of the target. Due to the difference in velocities 
between neutrons and y-rays the corresponding signals triggered 
the TAC at different times. Consequently the output signals of the 
TAC were essentiallysorted out into two groups separated by the 
time difference between the arrival of neutrons and y-rays to the 
Ge(Li) crystal. The scatterer-to-detector distance was sufficient 
to yield a good time separation between neutrons and y-rays in the 
TAC time spectrum. To reduce background y-rays a single-channel 
analyser (SCA) was used to set a window on the y-peak of the time 
spectrum. The SCA output signals were then used to gate the energy 
signals; both were fed to an amplitude-to-digital converter (ADC) 
running in a coincidence mode and interfaced to a Honeywell 316 


on-line computer. 
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The neutron flux was monitored by a fixed position Ne 213 
neutron detector employing time-of-flight techniques. The neutron 
Spectrometer has previously been described in detail (E1 69). 

Gamma ray spectra from the °®Fe(n,n'y) reaction were obtained 
at maximum neutron energies ranging from 0.89 to 1.92 MeV. The 
corresponding energies of protons falling on the Li target range 
from 2.60 to 3.60 MeV. "Scatterer-in" and "scatterer-out" runs 
were made at each energy for a predetermined charge collected on 
target. Peak areas of the 847 keV y-ray were extracted from the 
gated y-ray spectra. At each energy the intensity of the ground 
state neutron group of the ’Li(p,n)’Be reaction was also extracted 
from the neutron monitor spectrum. 

The efficiency of the Ge(Li) detector was determined using 
calibrated y-ray sources (Ro 72), and the neutron detection effic- 


iency was calculated using a Monte Carlo routine (see Appendix A). 


3.2 Data Reduction 


The quantity measured experimentally is the yield Yee sy) 
for a particular y-ray at a proton energy Ep and observation angle 
py with respect to the proton beam. It is desirable to describe 
this yield in terms of a y-ray differential production cross section 
at an effective neutron energy. Because of the nature of the experi- 
mental arrangement, this description must explicitly take into 


account the effect due to averaging the neutron energy and flux 


variation over the scatterer geometry. 
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Neutrons produced by the ’Li(p,n)’Be reaction are considered 
to come from a point source placed at the centre of a circular disc 
of radius R and thickness T. The y-yield for a particular Y-ray 


is given by: 


ed vea Sean 


| do, 
EpVy) = Ne ey | [1q(E,s0q)] Lge (Eqsey)] + Lexp(-u,t, 


Vv 
where No is the number of scatterer nuclei per unit volume in the 
sample, e, is the total efficiency of the y-detector including the 


solid angle, 1 CE 6) is the neutron intensity (cm* s') in the 


p 
volume element dV at a proton energy E, and neutron emission angle 
do 
6, with respect to the proton beam, pay Ene is the differential 


production cross section for the y-ray of interest at a neutron 
energy E, and angle 0, with respect to the neutron direction, Hy is 


the y-absorption coefficient and 2 is the distance in the sample 


nif 
through which the y-ray has to pass. The integration extends over 
"the active volume" of the scatterer which will be defined later. 


The first step is to calculate EAE »6,). This must take into 


p 
account the angular distribution and attenuation of neutrons in the 
scatterer. The ’Li(p,n)’Be reaction yields two neutron groups, no 
and ni, leading to the ground and first excited states of ’Be, 
respectively. For simplicity we will neglect in the following 


derivation the small contribution to the neutron intensity from the 


ny group. The differential cross section for the ’Li(p,no)’Be 
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reaction can be described in terms of a Legendre polynomial series 


(see Appendix A): 


Blep abs ccOsocny [3.2] 


If the neutron source is considered to be at the origin of a 
cylindrical coordinate system with the proton beam in the z-direction 


(fig. 3.2a), then I, is given by: 


Nel ) A.(E_) P; (cos anal ce 2 fouou 


where N is the number of ’Li nuclei per cm’, No is the number 


a) 
of protons falling on this target per second and oS aes tan™ (o/Z). 
Due to the small size of the sample, neutron attenuation is expected 
to be small. However, it can be accounted for in the determination 


of the factor N iN as described below. 


Lip 

Attenuation of y-rays in the scatterer must be calculated 
explicitly because typical y-ray absorption lengths are comparable 
to the scatterer dimensions. The expression for the y-ray path 
through the scatterer, ds can be determined with the aid of fig. 
3.2b. In the present experimental arrangement, where Wy = 90° and 
the scatterer-to-detector distance is large compared to the sample 


dimensions, the y-rays emitted within the solid angle of the detector 


are essentially perpendicular to the z-axis. Thus 2, is independent 
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FaGuUirene 3% 2 


(a) 


(b) 


Schematic drawing of a cylindrical scatterer 


and a point source located at S. 


Front-view of the scatterer with the z-axis 
perpendicular to the plane of the page and 


passing through the neutron source. 


af 


of z and is given by 


r| ~ 


& (os) = [R°- p* sin’ 6] - 9 cos ¢ [3.4] 
Eq. (3.1) then becomes 
do 
Y(Epo,) = ¢ | L1,(A,s042)]* Cat (E20) I+ Lexp(-u,2, (0.6)) 3 


e do dd dz i325 


o7 + 2 


Since both EH and ey are dependent on the value of o/z, 
it is clear from eq. (3.5) that the measured y-yield represents 
an average over both the neutron energy spread and the y-ray 
angular distribution. It is then possible to define this average 


production cross section as 


COPE nt = LY (E, )1 (| [1 (A;3042)ILexp(-u,2, (046))] 


o do do dz} [3.6] 
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This average differential cross section is associated with an 


"effective" neutron energy E. defined by 
Ey =f (C1) (Aj30.2)] + Lo, 4+ (E,)] + CE, (0.2)]0 do dz} 
x ({ Ot, (A,39.z)] loathe pedeedz yi fad 


where om nave) is a function that describes the dependence of 


the inelastic cross section on energy. It is taken to be (Bl 52) 


Onn (Ey) = (En-E ep) [3.8] 


The limits of integration over the po and z variables can be 
determined in the following manner. For a proton energy EO let 


th 


Bt be the neutron energy at 0° and Ee be the neutron threshold 


energy for the level decaying by the particular y-ray of interest. 


If a is the angle that corresponds to oa 


» then only neutrons 
emerging in the cone a will have enough energy to produce that 
Yv-ray. This cone will define what we called "the active volume" 
of the scattering sample. The angle a can be calculated for a 
particular level at each proton energy from the kinematics of the 
7Li(p3no)7’Be and the (n,n') reactions. A maximum value of 1/2 is 
used for a since backward-angle neutrons will not interact with 


the scatterer. Referring to fig. 3.3, the limits of integration 


are: 
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Side-view of the scatterer indicating the limits of 


integration over o and z variables. 
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The integration over the z-variable in eq. (3.6) can be 
carried out analytically (see Appendix B), whereas the integra- 
tions over p and ¢ have to be done numerically. The double 
integral of eq. (3.7) can be reduced to a single one over the 


angle 8, by using spherical coordinates. Thus 


| 


nea )raciol [+ f(0,) E,(0,) 46,1 xL fn (0, ) dey [3.9] 
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with the weighting function F(8,) given by 


atreene 3: 10} 


isc lees on (E, +01) - Rhee) es iney > a 
where & (0, ) is the distance travelled by a neutron emitted at 
an angle 0, The function F(6,) represents the number of inter- 
actions occurring at neutron energy EL 

Since the differential cross section given by eq. (3.6) 
“averages over the y-ray angular distribution", the total cross 


sections can be obtained from the relation: 
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A computer program was written to manipulate the raw data 
according to the above formulation (E172). The 
average differential production cross sections and effective neutron 
energies were calculated using eqs. (3.6) and (3.8), respectively. 
Legendre polynomial expansion coefficients for the ’Li(p,no)’Be 
reaction were obtained from a separate experiment (Appendix A). A 
correction for the small contribution from neutrons leading to the 
first excited state of ’Be was included in the determination of 
neutron intensity. The product Nia No was obtained for each run 
from the yield of the no group observed in the neutron monitor 
detector and the known ’Li(p,n,)’Be cross section. The neutron 
yield was corrected for detection efficiency and attenuation in the 
sample. The ratio of the No yield in two consecutive runs with and 
without the scatterer provided the attenuation factor. 

Experimental errors for the average differential cross sections 


evolve from uncertainties in the determination of the following 


quantities: 

) peak areas of y-rays and neutrons 

(ii) detection efficiency for y-rays and neutrons 

(iii) cross sections for the 7Li(p,n)’Be reaction 

(iv) number of scattering nuclei 
) y-ray absorption coefficients and neutron attenuation 
) integration over scatterer geometry. 


The overall errors in absolute differential cross sections were 
estimated to be < 25%. 
For the total cross section obtained from eq. (3.11), an 


additional error arises from the assumption that the present 
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geometrical arrangement averages over the y-ray angular distribu- 
tion. In order to estimate this error, let us write the differ- 
ential cross section as: 


do 
a (E28) = Ay [1+ a2 po (cos 8) [22123 


The above assumption implies that the integration of the second 
term of this equation over the scatterer geometry is negligible. 
To show that this is true the following integration was carried 


Out: 


I = | f(e,) Do (cos 6) d0,, J+ Ce, de]-" cong ec 


where F(6.) is the corresponding F(6, ) as given by eq. (3.9). 

The limits of integration were set according to the range of 6, 

as determined by the angle a of the neutron cone. It was found that 
the maximum value of I = -0.17 for sis = 100 keV above threshold 
and rapidly decreases to 0.01 with increasing energy. Assuming 

a value of 1.0 for a, the total cross section might be subjected 


to a maximum systematic error of 17% just above threshold. 


3.3 Experimental Results 


The timing characteristics of the y-ray spectrometer are 
indicated in the time-spectrum shown in fig. 3.4. The effect 


of the electronic time window on the observed y-ray spectrum is 
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Figure 


3.4 


Time spectrum from the TAC associated with the 


Ge(Li) detector (see fig. 3.1). 


00zI 


001 


000l 


YJgWNN ISNNVHD 


006 008 | 002 009 00S 5 
0001 
za 
c 
axg 
5 0002 
GH 
OO0€ 
006249 WOF= HIVd LHONS 
A®W S'Lexow 3 LY WNdld3adS 3WIL 1139 
(A,u’u) oe 000” 


TANNVHD Y3d SLNNOD 


34 


Shown in fig. 3.5. The suppression of events due to time-uncorrelated 
room background and neutron induced events in the detector is quite 
effective. The 431 and 478 keV y-rays are due to the ’Li(p,ny) and 
"Li(psp'y) reactions, respectively; the 847 keV y-ray corresponds to 
the well-known decay of the first excited state of °°Fe. 
The excitation curve of the 847 keV y-ray from °°Fe is shown 

in fig. 3.6, together with previous results obtained by conventional 
methods. It is quite clear that the present source-scatterer geometry 


performs a good energy average and yields reliable cross sections. 
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Figure 3.5 


Fe(n,n'y) spectrum obtained in a 51 min. run. The 
gated and ungated spectra were collected simul ta- 
neously. Peaks are labelled by their corresponding 
y-ray energies in keV and the material from which 


they originate. 
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Figure 3.6 


Excitation curve for the 847 keV y-ray from °°Fe. 
The average total cross section is plotted against 
the effective neutron energy. Error bars represent 
absolute errors. Previous results of the conven- 


tional technique are also shown. 
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CHAPTER IV 


APPLICATIONS: STUDY OF LEVELS IN**Sm AND 1°°Gd NUCLEI 


4.1 Experiment and Data Reduction 


The experimental arrangement described in § 3.1 was used. 
Scattering materials were obtained from Oak Ridge National 
Laboratory in the form of Sm,0, and Gd, 0,3 powders enriched to 
99.54 and 95.95% in !°*Sm and 1°°Gd, respectively. These powders 
have been compressed into a solid disc inside a nylon ring. 

Table 4.1 gives the weight and dimensions of both scatterers. 

Gamma-ray spectra following neutron inelastic scattering 
from }°*Sm and 1°°Gd were obtained at maximum neutron energies 
from 1.10 to 2.32 MeV. At each energy three separate runs were 
made: (i) a run with the scatterer (ii) a background run with 
an empty nylon ring, and (iii) a run with an Fe sample. The 
Fe run served as a check on cross section normalization. 
"Scatterer-out" runs were used to identify background y-rays 
observed in the time-gated spectra due to contaminant reactions 
created in the immediate surroundings of the scatterer. After 
normalizing "scatterer-out" runs and subtracting them from 
"Scatterer-in" runs, spectra due entirely to the scattering material 
were obtained. A typical spectrum obtained with the Gd scatter at 

m 


Eas = 1.71 MeV is shown in fig. 4.1. These runs were taken with a 


wider time gate admitting some of the y-rays generated in the Ge(Li) 
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TABLE 4.1] 


Weight and dimensions of scatterers 


Material Diameter Thickness Weight 
(cm) (cm) (g) 
evs 1.58 1.63 6.96 


(99.54% enriched in !°*Sm) 


Gd20 3 


1458 0.56 3.46 
(95.95% enriched in !°°Gd) 
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Figure 4.1] 


(a) 


160Gd(n,n'y) spectrum obtained in a 3 hour 
run. Contaminant peaks are labelled by their 
energies in keV and, whenever possible, the 


material from which they originate. 


*6°Gd(n,n'y) spectrum after subtracting a back- 
ground run collected for the same neutron flux. 
The peaks which are labelled by their energies 


in keV correspond to transitions in ?°°Gd. 
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detector. In fig. 4.la, the 596 and 690 keV y-rays are due to 
neutron interactions with the Ge crystal while the 847 keV y-ray 
is due to neutron interaction with Fe materials around the target. 
None of these lines appear in the subtracted spectrum (fig. 4.1b). 
Gamma-rays for which energies have been assigned in fig. 4.1b are 


attributed to transitions in +©°Gd. A similar spectrum for !°*Sm 


40 


is shown in fig. 4.2. Due to the very high background below 500 keV 


from the ’Li(p,p'y) and ’Li(p,ny) reactions, (see fig. 3.5) no 


quantitative analysis was carried out in this region. 


Peak centroids and areas of the observed y-rays were extracted 


at each energy using a peak fitting program (Ea 72). Gamma ray 
energies were obtained from peak centroids and a calibration line 
that was independently determined using radioactive sources of 
known y-ray energies. Production cross sections and effective 
neutron energies were calculated employing the method of analysis 


given in § 3.2. The y-ray absorption coefficients for Sm203 and 


Gd.0, were estimated using the tables of Davison and Evans (Da 52). 


Branching ratios for y-rays associated with a particular level 
have been determined from the production cross sections of these 
y-rays at all energies. The inelastic cross section for a given 
level was obtained by summing the cross sections for all y-rays 
originating at this level and subtracting from this value the 


cross sections for all cascades feeding this level. 
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Figure 4.2 


154Sm(n,n'y) spectrum after subtracting a back- 
ground run collected for the same neutron flux. 
The peaks which are labelled by their energies 


in keV correspond to transitions in ?°*Sm, 
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4.2 Statistical Theory Calculations 


As mentioned in Chapter II, the cross sections predicted by 
the statistical theory represent an average over many levels in 
the compound nucleus. Since nuclei of interest in this work 
(*°*Sm and 1®°Gd) have a high mass number, the excitation energy 
of the CN is then great enough (~ 8 MeV) so that the density 
of levels is large. In addition, the present experimental 
arrangement utilizes a neutron beam of a wide energy spread, 
thus resulting in the excitation of many CN levels. Accordingly, 
the measured inelastic cross sections are energy averages and 
therefore the statistical theory can safely be applied to predict 
such cross sections. 

The nuclei '°*Sm and *©°Gd fall in the deformed region of 
nuclei. Thus one would not expect to be able to fit experimental 
data assuming a spherical optical potential. However, the extensive 
work of Auerbach and Moore (Au 64) has revealed that one can 
reproduce, quite well, both elastic and inelastic neutron scatter- 
ing data on deformed nuclei in the region between Ta and U using 
a spherical, local, spin-dependent optical potential. Similar 
results for Ta, Re and Pt were reported by Smith et al. (Sm 68). 

A computer code has been written (Sh 72) to calculate neutron 
transmission coefficients. We have employed a spherical local 


optical potential of the form 
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: bide eae 
U(r), = - V TUR Oe say) + 4i a; W ap Efdrsr,s 2a.) + 


( sae V 1d rey Ye Gael eae [4.1] 
mc s yr dr Py OST TS . 


where V, W and ee are the depths of the real, imaginary, and spin- 
orbit wells, respectively, and m_. is the mass of the pion. 
F(rsr 2a) are taken to have the Saxon-Woods shape 


1 1 


f(r.r,,a) = {]+exp ieee sae 


where ee and a are the radius and diffuseness parameters respectively. 

In what follows, the depths are given in MeV and the geometrical 

dimensions in fm. 

Four "global" sets of potential parameters are available from 
analysis of neutron data over a substantial energy range and for a 
series of nuclei. These are: 

(a) Perey and Buck non-local potential (Pe 62) that fits neutron 
scattering data from 1 to 25 MeV. An equivalent over-al] 
local potential has been obtained by Wilmore and Hodgsen 
(Wi 64); this is specified by the parameters 
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Where E is the incident neutron lab energy in MeV. 

Moldauer potential (Mo 63): The Perey-Buck potential is not 
able to fit the neutron S-wave strength function around A = 100. 
Moldauer (Mo 63) was able to show that this could be corrected 
by allowing the absorbing part of the potential to peak sharply 
just outside the nuclear surface. Analysis of available data 
for neutrons of energy below 1 MeV led to the potential (given 


according to the geometry specified in eq. (4.1)) 


V = 46.0 hae p25 age ts 0.62 
W = 14.0 a 1.38 a. = 0.24 
Ve = 0.0 


Rosen potential (Ro 66): This potential was determined by 
fitting 14 MeV neutron elastic-scattering data on a wide 


range of nuclei. It is specified by the following parameters 


V = 49 - 0.33 E ena 1.25 ani as 0.65 
Wo= 5.75 a aa TeZo ae 5 Ory 
ve =. 5,5 OO Veco Qe 0.65 


The energy dependence of the real part of the potential was 
obtained from analysis of proton polarization data at differ- 
ent energies. 

Becchetti and Greenlees potential (Be 69) which was determined 


from analysis of neutron data for energies <14.5 MeV and nuclei 
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with A > 40. They utilized differential cross sections, 
polarizations, reaction cross sections and total cross sections 
and obtained excellent over-all fits. The optimum neutron 
parameters were found to be (for energies < 7.1 MeV) 


N-Z 


Vv = 56.3-0.32E - 24 SA Fea aimlel? a. = 0.75 

Woe 13 -0.25E- 12 Se <Re 7: Fee Diy 
; A 01 ; j 

Ne = 6.2 eS a, = 0.75 


where N and Z are the number of neutrons and protons in the 

nucleus, respectively. 

It should be noted that all these potentials yielded a reasonably 
good overall fit to neutron elastic scattering data down to ~1 MeV 
incident energy, provided that one takes account of compound elastic 
Scattering. 

In order to test the sensitivity of inelastic cross sections to 
these different potential sets, we have compared the measured and 
calculated cross sections forthe 847 keV state of °°Fe (fig. 4.3). 
The experimental results are in good agreement with the curves 
predicted by the potentials (a), (c) and (d) but in poor agreement 
with that predicted by the potential of Moldauer. Thus, it seems 
that Moldauer's average potential does not describe iron very well. 
Similar results were reported by Gilboy and Towle (Gi 65). 

Since the potentials (a,c and d) give similar results, we 
have arbitrarily chosen the Rosen potential for calculations carried 


out for Sm and Gd nuclei. Inelastic calculations reported in this 
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Figure 4.3 


Inelastic scattering cross sections of neutrons 
leading to the 847 keV state in °°Fe. Theoretical 
curves were calculated using the potentials of 
Wilmore and Hodgson (Wi 64), Rosen et al. (Ro 66), 
Becchetti and Greenlees (Be 69) and Moldauer 


(Mo 63). 
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thesis were carried out using Moldauer's code "NEARREX" (Mo 64a). 
The width fluctuation correction was applied and the (n,y) exit 
channel was taken into account. Input parameters required for 
this channel were obtained from the tables of Goldberg et al. 

(Go 66) and Gilbert and Cameron (Gi 65a). However, the effect 


of including the radiative capture channel was negligible. 


4.3 The Sm Results 


Level positions and decay modes are summarized in table 4.2. 
These were obtained from an error-weighted average of the results 
at all energies. The values obtained at the different energies 
showed good agreement in each case. The levels will be discussed 
in order of increasing energies. 
4.3.1 The 921, 1202 and 1286 keV levels 

The 921 keV level has been assigned J" = 1” by Veje et al. 
(Ve 68) from a study of (d,d') excitation cross section systematics. 
The observed decay modes (table 4.2) are consistent with this 
assignment. The present (n,n'y) excitation function is shown in 
fig. 4.4 together with the excitation cross section for a1 level 
predicted by the statistical compound nuclear theory. It was not 
possible to fit the observed yield more than 200 keV above threshold. 
As the two y-rays assigned to the 921 keV level have the same 
excitation curve shape, and as the present branching ratio determin- 
ation is in good agreement with that of D'Auria et al. (Da 71) from 


B-decay studies, it is most likely that the cross section 
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Excitation energies and branching ratios 
of levels in ?°*Sm 


Se ance Lee ee eee Pe ee eee eee eee ees 
Initial state Final state Ey Branching ratios 
(keV) (keV) (keV) Present work D'Auria et al. 


921 4.1 g.S. 921 42 +1 40 
82 839 58 +1 60 
1012-4 1 82 930 62a 1 58 
267 745 38 +1 42 
1099 + 1 82 1017 100 100 
1178 + 1 qis% 1178 DG? 28 
82 1096 40 + 2 40 
267 911 34 + 2 32 
1182 + 2 267 915 100 
12027201 82 1120 (50) 
921 281 (50) 
1286 + 1 82 1204 (75 
921 365 (25) 
1338 + 1 82 1256 46 + 3 
267 1071 54 + 3 
1440 + 1 dys. 1440 ho + 2 44 
82 1358 522 51 
267 1173 Br 2 5 
1475 +1 82 1393 90 
1099 375* 10 
1515 +1 82 1433 100 
1539 + 1 82 1475 oo en 72 
267 1272 D5 4 2 14 
1012 527 14 
1585 + 2 82 1503 7 a2 33 
267 1318 64 + 2 67 
921 664 29 + 2 0 
1755 + 2 921 834 ~ 75t 37 
1012 743 ~ 25 30 
1585 170* 33 


Intensity ratio (see text). 
* These y-rays could not be observed in the present experiment 
(see text). 
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Figure 4.4 


Excitation curves for the 921 keV level and the 
1120 and 1204 keV y-rays which depopulate the 1202 
and 1286 keV levels, respectively. The probable 
contribution from the latter two levels to the 

921 keV level, as discussed in 84.3.1, is also in 
indicated. The experimental points for the 92] 

keV level above a neutron energy of 1.5 MeV contain 
a ~ 15% correction due to feeding from the 1585 and 


1755 keV levels. 
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discrepancy is due to unaccounted-for cascade feeding from a level 
or levels between 200 and 500 keV higher in excitation. The most 
likely candidates are the 1202 and 1286 keV levels, as discussed 
below. 

The presence of a level at 1202 keV excitation is indicated by 
the excitation function of the 1120 keV y-ray (fig. 4.4). This is 
most likely the J" = 0° level observed in the (t,p) reaction at 
1218 keV. Also the excitation function of the 1204 keV y-ray 
shown in fig. 4.4 indicates an origin at 1286 keV excitation. This 
is most probably the level observed via the (p,p') reaction (Ke 64) 
at 1295 keV and at 1299 keV via the (t,p) reaction (Bj 66). 

It should be mentioned that at the early stages of the present 
analysis it was thought that the 1120 and 1204 keV y-rays belonged 
to a level at 1204 keV which is quoted inthe paper of D'Auria et al. 
(Da 71). However, as the two y-rays appear to have different 
threshold energies (fig. 4.4) it has been concluded that the two 
y-rays belong to separate levels as mentioned above. 

The 1202 and 1286 keV levels can decay to the 921 keV level 
via 281 and 365 keV y-rays, respectively. Both y-rays have been 
observed in the previous y-decay work (Da 71). The 281 keV 
transition was assigned to two other levels, but the assigned 
transition strengths in the decay scheme only account for some 
75% of the total observed strength; the 365 keV y-ray was not 
assigned to any level. A 365 keV y-ray was observed in the present 


work with an approximate threshold of 1.3 MeV but no quantitative 
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yield measurements could be obtained ( 84.1). 

The 1120 keV y-ray excitation curve is in agreement with the 
predicted shape for a ae level, but is only approximately half 
the expected strength. The 1204 keV y-ray yield can be fitted by 
either a J =3 curve scaled by a factor of 0.75, or J=4 scaled by 
1.5. The observed branch to the 921 keV 1 level and the excita- 
tion of the level via the (t,p) reaction then strongly suggest an 
assignment of J" = 3° for the 1286 keV level. The expected 
contribution to the 921 keV level from these two states then brings 
the experimental yield curve into fair agreement with the expected 
excitation function. 
4.3.2 The 1012 and 1099 keV levels 

The 1012 keV level has been observed in several reactions 
(Ke 64, Bj 66, Se 66, Ve 68, Da 71), and has been assigned J's 
3 (Ha 63, Ve 68). The present experimental yield curve for this 
level is shown in fig. 4.5a. Also shown are various yield curves 
which also include the expected contribution from the 1372 keV 
a level and the known (Da 71) decay scheme for the 1372 keV level. 
The experimental points indicate a preferred assignment of ae 
in good agreement with the previous results. 

The 1099 keV level decays entirely to the 82 keV (2°) level 


“ (Bj 66, Ve 68). The experimental 


and has been assigned to 
yield curve of the 1099 keV level has been corrected for cascade 
feeding from the 1475 keV level using the measured 1475 > 82 keV 


transition yield and the known (Da 71) decay scheme for the 1475 keV 
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Figure 4.5 


Excitation curves for the 1012 (a) and 1099 (b) keV 
levels. The calculated curves for the 1012 keV 
level also contain the expected contribution from 
the 1372 keV ae level, and the experimental points 
above 1.7 MeV have been corrected for feeding from 


the 1707 and 1755 keV levels. 
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level. The resulting yield curve (fig. 4.5b) indicates a preferred 
assignment of te Or: in agreement with the previous results. 
4.3.3 The 1178 and 1182 keV levels 

A close doublet has been observed at approximately 1.18 MeV 
excitation with suggested spins of ao and 5 (Yo 65, Se 66, Ve 68). 
The present experiment has clearly resolved the members of this 
doublet to be at 1187 and 1182 keV. The yield curves (fig. 4.6) 
indicate a 2° and 5” assignment for the 1178 and 1182 keV levels, 
respectively. D'Auria et al. (Da 71) have also been able to 
measure the decay modes of these two levels. 

The 636 keV y-ray, assigned by Seaman et al. (Se 66) and 
D'Auria et al. (Da 71) to the decay of the 1182 keV level, was 
not observed in the present work below a neutron bombarding 
energy of about 1.8 MeV and so cannot belong to the 1182 keV level. 
This y-ray could be due to a transition from a level at 1818 keV 
to the 1182 keV level; a level at 1811 keV has been reported in the 
(d,d') work (Ve 68). 
4.3.4 The 1338, 1440 and 1707 keV levels 

The excitation function of the 1071 and 1256 keV y-rays (fig. 
4.7a) indicates a common origin from a level at 1338 keV excitation. 
This probably corresponds to the level observed by Kennefick and 
Sheline at 1344 + 10 keV. The observed decay scheme and excitation 
curve indicate an assignment of q* for the 1338 keV level. 


The 1440 keV ae level decays via 1440, 1358 and 1173 keV y-rays 
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Figure 4.6 Excitation curves for the 1178 and 1182 keV levels. 
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Figure 4.7 


Excitation curves for the 1338 (a), 1440 and 1707 
(b) keV levels. The square points in the excitation 
curve of the 1440 keV level have been corrected for 
the contribution from the 1707 keV level using the 
previously measured (Da 71) branching ratio and the 


+ j : : 
expected 4 excitation cross section. 
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to the ground, 82 keV and 267 keV levels, respectively. The 1440 keV 
y-ray has also been assigned to the decay of the 1707 keV level 
which has a spin of either 3. or 4* (Da 71). The contribution to 
the yield curve due to the 1707 keV level was estimated using its 
known decay scheme (Da 71) and the expected cross section for both 
3” and 4". However, it was possible to obtain a consistent 
branching ratio determination for the 1440 keV level only if the 
spin of the 1707 keV level was assumed to be 4”. The yield curve for 
the 1440 keV level is then in reasonable agreement with the predicted 
2* excitation function (fig. 457 by 
4.3.5 The 1475, 1662 and 1515 keV levels 

The 1393 keV y-ray belongs to the decay of both the 1475 
and 1662 keV levels (Da 71) which have been tentatively assigned 
emeel ave 6S) cand J". = 4* (Ve 68, Yo 65) respectively. The 
assumption of two origins for the 1393 keV y-ray is supported by 
an observed 2 keV energy shift between bombarding energies of 
EY = 1.6 and 2.1 MeV in the present work. 

The yield curve of the 1393 keV y-ray shown in fig. 4.8a is 
compared with the expected contribution from the two levels using 
the above assignments. Reasonable agreement is obtained for both 
the 1475 and 1662 keV levels using the expected yield curves for 
1” and 4” levels, respectively. 

At neutron energies above 1.6 MeV,a 1433 y-ray was observed 


(fig. 4.8b). The yield curve indicates that this y-ray could 
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Figure 4.8 


(a) 


(b) 


Excitation curve of the 1393 keV y-ray 
depopulating the 1475 ana 1662 keV levels. 
The predicted curves for the 1475 keV levels 
have been adjusted to account for the 
(unmeasured) 10% decay to the 1099 keV 
level. The y-ray from the 1662 keV level to 
the 82 keV level was weakly observed in the 
work of D'Auria et al. (Da 71). 


Excitation curve of the 1515 keV level. 
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belong to a level at either 1433 or 1515 keV excitation. No 
previously reported level would correspond to an excitation of 
1433 keV, but a level, tentatively assigned J = an has been 
observed at 1522 keV by Kenefick and Sheline (Ke 64), and so 
the 1433 keV y-ray is assumed to belong to that level. No other 
y-ray could be assigned to the level. From a comparison of the 
observed yield with that expected for various possible spins, 
the 1515 keV level can have a spin of 1 or 3. Of these, the 
decay mode is more consistent with J = 3. 
4.3.6 The 1539 and 1585 keV levels 

Transitions to the 82, 267 and 1012 keV levels have been 
observed from a level at 1539 keV. Although D'Auria et al. (Da 71) 
have proposed this level as a new one, it is more likely to be the 
same level reported in the (p.p') work (Ke 64) at 1547 keV (note 
also the 7 keV difference in the case of the 1515 keV level). In 
the present work the transition strength to the 1012 keV level was 
not measured and consequently the yield curve shown in fig 4.9a 
has been corrected for this decay using the known branching ratio 
(Da 71). The measured excitation function for the 1539 keV level 
can be fitted with either J = 1 or 3; the decay scheme then 
indicates a preferred assignment of 3. 

The level at 1585 keV has been assigned J" = 3° by Veje et 
al. (Ve 68). The present results (fig. 4.9b) are in good agreement 


with that assignment. The 664.6 keV y-ray was assigned by D'Auria 
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Figure 4.9 


Excitation curve of (a) the 1539 and (b) the 
1585 keV levels. Branching ratios of y-rays 
from the 1539 keV level are those of D'Auria 


et al. (Da 71). 
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et al. to the decay of '°°Sm, produced by neutron capture. However 
in the present work the excitation function of the 664 y-ray was 
entirely consistent with the decay of a level in '°*Sm at approx- 
imately 1.58 MeV. Consequently it has been assigned to the decay 
of the 1585 keV level. 

4.3.7 Other levels 

The 743 and 834 keV y-rays from the 1755 keV level were 
observed at neutron energies greater than 1.8 MeV. There was no 
evidence at any energy for the existence of a 170 keV y-ray assigned 
to the 1755 keV level by D'Auria et al. In the latter work, y-rays 
within 1 keV of this energy were also assigned to decays in 1°!Nd 
and +°°Sm, which may account for the discrepancy with the present 
work. 

The close proximity of the 743 keV y-ray to the much stronger 
745 keV line made it difficult tc obtain reliable estimates of the 
743 keV y-ray yield. However, the determined branching ratio 
obtained in the present work is in fair agreement with the intensity 
ratios obtained by D'Auria et al. The measured total yield curve 
(fig. 4.10a) suggests an assignment of = for the1755 keV level. 
This is consistent with the observed decays to the 1012 keV 3 level 
and the 921 keV 1 level. 

The 753 keV y-ray was observed for neutren energies above 
1.8 MeV. This y-ray does not correspond to a decay from any 
previously observed level, and most likely indicates the existence 


of a new level at an excitation of either 1674 or 1765 keV. The 
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Figure 4.10 


Excitation curve of (a) the 1755 keV level and 
(b) the 753 keV y-ray. The calculated yield 
curve for the 753 keV level indicates the two 


possible level origins. 
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excitation function then suggests a spin of 4* or 0° if the level is 
at 1674 keV or 3 if it is at 1765 keV.(fig.4.10b). The decay mode 
is most consistent with an assignment of the 753 keV y-ray to a 
level at 1674 keV with a probable J" = 0°. 

A 1550 keV y-ray was observed in the present work only at 
neutron energies greater than 1.8 MeV. Most likely this transition 
and the 636 keV (see §4.3.3) belong to a level at 1817 keV. If this 
level corresponds to that observed by Veje et al. (Ve 68) at 1811 
keV with J™ = 3 or 4*, then the present decay mode is consistent 
with J = 4". D'Auria et al. (Da 71) has assigned the 1549 keV 
y-ray to the decay of a new level at 1549 keV. 

4.3.8 Level scheme and band structure of !°*Sm 

The level decay scheme of 1°*Sm up to an excitation of 1755 keV 
is shown in fig. 4.11. Rotational bands built on octupole, B- and y- 
vibrational modes were proposed in +°*Sm by many investigators. 
Members of these bands are also shown in fig. 4.11. We have been able 
to clearly resolve the third member of the octupole band at 1182 keV 
from the second member of the 8-band at 1178 keV and to determine 
their spins to be 5 and oa respectively. The level at 1539 keV was 
identified as the 3 member of the y-band. This is consistent with 
both the energy systematics of this band (Ku 67) and the observed 
decay mode and inelastic cross sections. The third member of the 
@-band shown at 1372 keV excitation was not included in the present 
study as the predominant decay mode is too low in energy to obtain 
reliable information. 


Most of the unclassified levels shown in fig. 4.11 were reported 
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Figure 4.1] 


Level scheme of '°*Sm excited states up to 1755 
keV; known members of various bands are also 
indicated. The excitation energy of the 8* level 
shown at 906 keV is based on the energy of the 

8° > 6° transition measured by Diamond et al. 

(Di 72) and the position of the 6° level reported 
by D'Auria et al. (Da 71). The decay mode of the 
1372, 1475, 1539 and 1707 keV levels are from 
D'Auria et al. (Da 71). 
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in the (p,p') work but were not observed in other experiments. It 
was possible to populate them in the present reaction study and 
observe their de-excitation y-rays. In addition, possible spin 
assignments for these levels were made and indicated in fig. 4.11. 
The 1 and 3. levels at 1475 and 1585 keV could possibly be members 
of a K = 1 octupole vibrational band. It must be anticipated that 
any additional classification will be very difficult. 

In general, the level spin assignments and branching ratio 
determinations obtained from the present study are in good agree- 


ment with previous information. 


4.4 The Gd Results 


The excitation energies and decay modes for levels of *°°Gd up 
to 1600 keV excitation are summarized in table 4.3. These are 
error-weighted averages of the results obtained at different 
bombarding energies. 

It seems appropriate, before going into detailed discussion 
of results, to mention that the calculations of inelastic cross 
sections were carried out using four sets of transmission 
coefficients which correspond to the four potential parameters 
listed in §4.2. This was done in order to check the sensitivity 
of inelastic calculations in the deformed mass region to different 
spherical optical model potentials. Predictions based on the 


potentials of Wilmore and Hodgson, Rosen et al., and Becchetti 
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TABLE 4.3 


Excitation energies and branching ratios 
of levels in *°°Gd 


a S 


Initial state Final state Ey ; e a 
rancnin ratios 
(keV) (keV) (keV) i 
989 + 1 gG.S. 989 46 + 1 
75 914 54 + ] 
1058 + 1 75 983 83 + ] 
249 809 7 
1071 +1 75 996 36 41 
249 822 64 + ] 
1149 + ] 75 1074 34+] 
249 900 66 + 1 
1225 + ] a.$3 1225 39 +1 
75 1150 61 +1 
1262 + 1 75 1187 66 + | 
249 1013 34+] 
1290 + 1] 75 1215 54 + Qt 
249 1041 46 +2 
1351 +1 Gs Se 1351 29 fty3 
75 1276 78 + 3 
1378 + 1 qc 1378 100 
1381 +2 Ges 1381 100 
1464 + 1] 75 1389 
249 1215 
1569: tn2 O53 1569 
75 1494 
989 | 580 
1587 +2 Gabe 1587 
75 1512 
1599 +2 Gass 1599 
75 1524 


t These are the values determined at a < 1.4 MeV (see text). 
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and Greenlees (§ 4.2) were similar to within 10%. The "reasonable" 
agreement obtained between measured and calculated cross sections 
for the Sm (§4.4) and Gd (shown below) nuclei indicates that any one 
of these three potentials is good enough for inelastic calculations, 
even in the case of highly deformed nuclei. The Moldauer potential 
has yielded good fits to excitation functions of positive-parity 
States but underestimates, by a factor of ~2, the cross sections for 
levels of negative parity. In fact, the use of Moldauer transmission 
coefficients has improved the fits, especially just above threshold 
(in the case of positive parity states). Moldauer potential was 
designed to fit the s-wave neutron strength function and elastic 
Scattering data for energies below 1 MeV. This probably explains 
why it yielded better fits at low energies for positive parity 
SLaves hig.< 4.12 displays tne transmission functions calculated 
using both Rosen and Moldauer potentials. The difference between 
the values of To» especially for odd partial waves, is quite 
obvious. 

In the following we will discuss the levels of *®°Gd in order 
of increasing energy. 
4.4.1 The 989 keV level 

This level has been observed in both (p;,p') and (d,d') studies 
(Sh 66, Bl 67), and has been assigned Iss oe from cross section 
systematics of the (d,d') reaction (Bl 67). The observed decay 


mode (fig. 4.13) is consistent with either jg. 21> or 2). No 
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Figure 4.12 


The neutron transmission functions of '°®°Gd 
for various partial waves using Rosen and 


Moldauer potentials. 
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Figure 4.13 


Excitation curve of the 989 keV level. 
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branch to the 4° level of the ground state band could be identified 
in the present results. The strength of such a transition is 
expected to be very small (3%). The measured neutron inelastic 
cross sections are shown in fig. 4.13 together with theoretical 
predictions using Rosen and Moldauer transmission coefficients. A 
good fit is obtained for J" = a in agreement with previous 
assignments. This level has been identified by both Sheldon (Sh 66) 
and Bloch et al. (Bl 67) as the band head of the y-band. 
4.4.2 The 1058 keV level 

Sheldon (Sh 66) has reported a level at 1060 keV. If the 
989 keV level is the o* member of the y-band, then the ch member 
Should lie at ~1057 keV. On this basis, Sheldon has tentatively 
proposed a spin of 3° for this level. Since the (d,d') reaction 
preferentially excites 3. states (and not ate this level has not 
been observed by Bloch et al. (Bl 67). In the present work the 
y-rays of 983 and 809 keV have been associated with a level at 
1058 keV. Its decay mode and excitation curve (fig. 4.14) are 
consistent with an assignment of J = 3. From potential energy 
systematics, this level has been identified as the a member of 
the y-band. 
4.4.3 The 1071 and 1149 keV levels 

The level at 1071 was weakly excited by the (d,d') reaction 
(B1 67) and tentatively assigned J" = 2*. It has not been 


reported in the (p,p') work (Sh 66). Two transitions (996 and 
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Figure 4.14 


Excitation curve of the 1058 keV level. 


(mb) 


oe 


300 


200 


100 


1058 keV level 1058 


——— ROSEN POTENTIAL 
——_——=—— MOLDAUER POTENTIAL 


249 
IAS, 


83 


5 tee 


7 ot ‘ Se os Aa as. v 
- i ro - ee A 2 a ~ 
Tite ote a: ea One | 


. JAITTOY 
JAITHATOS 13 


822 keV) from this level to 2* and q* levels of the ground state 
band have been observed in the present study (fig. 4.15). If 

this level has a spin of an as suggested by Bloch et al. (Bl 67), 
it should decay by a relatively strong branch to the ground state. 
Assuming that the 1074 keV y-ray observed in the present spectra 
corresponds to such a pee oF transition (note that its energy is 
off by 3 keV), it was found that the magnitude of cross section 
for this level (combining the three y-rays) is only 40% of the 
theoretical values for a spin of Ae Moreover, the shape of the 
excitation curve indicates a high spin state. Therefore the 

1074 keV y-ray has been assigned to the decay of the 1149 keV 
level as discussed below. Consequently, the possibility of Jie 
2° for the 1071 keV level was ruled out. The measured inelastic 


cross sections were then compared with theoretical predictions for 
+ 


7] 


ied. Although the shape of the experimental excitation function 


is consistent with the one calculated using Rosen potential, only 
60% of the strength is observed. The use of Moldauer transmission 
coefficients yielded a good agreement with measured cross sections 
only for the region about 400 keV above threshold. If this level 
sea 4" then the transition to the 6’ member of the ground state 
band is expected to be very weak and was not observed. 

The level at 1149 keV has previously been identified as the 
4* member of the y-band (Sh 66, Bl 67). The 1074 and 900 keV 
y-rays observed in the present experiment have been assigned to 


the decay of this level to 2* and 4* members of the ground state 
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Figure 4.15 


Excitation curve of the 1071 keV level. Theo- 
retical predictions are for J” = q* using Rosen 


and Moldauer transmission coefficients. 
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band, respectively. Fig. 4.16 shows the decay mode and inelastic 
cross sections for this level. It is clear that the measured and 
calculated cross sections (for J" = 4*) show a situation similar 
to that encountered in the case of the 1071 keV level. It is 
also noticed that the branching ratios of y-rays from these two 
levels are the same. 

The 1149 keV level can decay to the 6" member of the ground 
State band by avery weak y-ray of 633 keV energy. Although such a 
y-ray was observed, we feel uncertain about assigning it to such 
a transition due to the relatively high background in this energy 
region. In addition, its excitation curve shows a sudden rise at 
neutron energy around 1.7 MeV which indicates that it might be 
a transition from a level at 1691 keV to the 1058 keV state; a 
3 level at 1688 keV has been observed in the (d,d') work. 
Regardless, the strength of this 633 keV y-ray below EY = 1.7 MeV 
is too small to account for the difference between measured and 
calculated excitation functions. 

4.4.4 The 1225 and 1290 keV levels 

The existence of these states is well established by previous 
studies.(Sh 66, Bl 67). A tentative assignment of 1 has been made 
by Bloch et al. (Bl 67) for the 1225 keV level. The present 
results (fig. 4.17) confirm this assignment. 

The 1290 keV state has been observed to decay to the 2* and 


4* states of the ground state band via 1215 and 1041 keV y-rays. 
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Figure 4.16 


Excitation curve of the 1149 keV level. 
Theoretical predictions are for a= 4X 
using Rosen and Moldauer transmission 


coefficients. 
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Figure 4.17 


Excitation curve of the 1225 keV level. 
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Branching ratios determined at different energies have indicated 
a slight increase in the strength of the 1215 keV y-ray at Ee 
1.5 MeV. This can be noticed in the excitation curve of the 
1215 keV y-ray shown in fig. 4.18. The 1215 keV was partially 
assigned to the decay of the previously reported 3 level at 
1464 keV (Bl 67). The observed decay mode and inelastic cross 
sections for the 1290 keV state (fig. 4.18) are consistent with 
a J" = 3° assignment, in agreement with previous results (Bl 67). 
4.4.5 Other levels 

The observation of y-rays which could not be fitted into the 
previously known level scheme of !©°Gd indicates the existence of 
new levels. Therefore, we propose levels at 1262, 1351, 1378, 
T1381, 1569, 1587 and 1591 ke¥- 

The 1262 keV level was observed to decay to the ae and 4* 
states of the g.s. band via 1187 and 1013 keV y-rays, respectively. 
Fig. 4.19 displays the excitation functions of these y-rays. The 
decay mode and inelastic cross sections suggest a spin of q* for 
the 1262 keV level. 

Two y-rays with energies of 1351 and 1276 keV were assigned 
to the decay of a level at 1351 keV. Fig. 4.20 shows the excitation 
curve of this level. A good agreement between measured and 
calculated cross sections was obtained for J" = 1. 

A y-ray doublet of 1303 and 1306 keV was observed at neutron 


energies around 1.4 MeV. Fig. 4.21 displays the excitation function 


76 


‘ cor _ } 
" ie See ite ae as 
-- BBSeIToNT oven geteiedS IgsistTto Je DBA Preiss 2dtte pation i 


Sy Fecuety Vor efSf edt to dd pnsyde .5H3 hi 26S a ont siete 8 


: bs - . ‘4 e 5 eee! ay 
' oi es + = #"% i} ‘24 & { F &3 : 7 9. 2 Ai réti prep i For od f 8D 2taT ne 2.f 
; S.-M TSIa ew Vad FST eit aT Uh ork ii “nwotte V6I-v Ved 21 PY 
—, 2 , a) = = 5 ae 7 
t ih lovsi £ bsivegsy vleuotvexd ont te Xe086 SH9 + Ot bang fe 4 
b WOT | 
; ! : : of : 
‘ ~~ * we : 7 : » \ he ‘ 3 . oe ren) 
ID QISESi Se) LAS SPOON NYSvon OovysecO SAt af eo a 
1 C tie, | fl 
Ait 4 ; Pega i < (RQ { iy aay i? Voy GeSk ot Yar 
=" . i Ps 
t \a bay 2s be f 2vior 1G ' wis ? aS fT: a dt ft¢ inp ee eh 
é fovey ToASD re 
dt ofai boettit ad ton bios dathw eveteu-t0 nohteyeedde Saey ASL 
f ipa of Ay 
Povommaterxs ord astearbnt bo°*' to onedbe TSvshiawom Yano rerg 
j C : : of Ls y : Z 
ra rat ener xf Lit. s) tao ~ ee et es " ne 
(hl . feet SAS! Fe efevs!l s2cdosg sy 2s iofsraity, -eteyal We ie f 
1 = ; , | 1 ru +" ve j 
; | es fee" bag Neel, vaet . fe set 
: \ 
ian 4 5 ee Pe apn sy oe ¥. CF @ eid 
i FY tins S edd 03 YAV30 OF  pOVIS edo Wa lS oy Yor ! cas i oF Z 
: f ee is M4 4: 
a a P* yy a a rine - | , 
« uf SY (229GeSF § . eva V Vox £10) brs - (8! i 5 ee bined, Sait e ‘silt aia zotese ¥ 
aay .2yerny seals Fo enetssnut notice rroyes sate ayetgarb QT. be 
4s ie 7 ; 7 i ; in { 
sa? 2B To nide 6 t2oppile Snotiss2 22049" visoatgnt bing. ahem oe 
. | i, ve | “Fave Vea ‘sosr oat 
i vale 


: f a ; ia : 
bongfe2b, vow Yan dy 37) DAS, fee to sl dis ayy Wr a 


“otdettoxs at dione 0S) int? cved) Fae Je. Viaush nto soa vb 7 
; ; ving _ 


bis bavuzbam nosy A ad tnambarg 16 ‘boog A ” [ovat atid % ow io 
hyo! is 
y} yo. Pre Oy vot beniatde Cian 20010. badietdten 
: | Agi P _ 


a i" 7) alive by bavicaus 25W Vad: 908 ne 088 4 bh wb yer A 
| ke cv 


=e ee 


7 fottonut foltstto ‘9 sad aveligetb ee rs OFT . 
ee 


1s 


M1 ¢ 


"© 1290 keV level 2. 
0) /® ADS) bev y-ray 


samme ROSEN POTENTIAL, 5c 
mmo MOLDAUER POTENTIAL 


77 


Figure 4.18 


Excitation curves of the 1215 keV y-ray and 
the 1290 keV level. 
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Figure 4.19 


Excitation curves of the 1187 and 1013 keV 
y-rays depopulating the 1262 keV level. Solid 


points are the total inelastic cross sections. 
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Figure 4.20 


Excitation curve of the 1351 keV level. The 
solid line is the theoretical prediction based 


on Rosen potential and for J" = 1°. 
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Figure 


4.2] 


Excitation curves of the 1303 and 1306 keV 


y-rays 
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of these y-rays. They were assigned to the decay of two levels 

at 1378 and 1381 keV to the 2° state of the g.s. band. The 
measured cross sections for the 1378 keV level can be fitted with 

a 1 curve scaled by a factor of 0.8 or 0° curve scaled by a 

factor of 1.4. The cross sections for the 1381 keV level could 

not be fitted with either 0” or 1 curve. No other y-rays could be 
assigned to the decay of these levels. It should be mentioned 

that because of the close proximity of the 1303 and 1306 keV y-rays 
the extracted yield of the weak 1306 keV y-ray was subjected to 
large statistical uncertainties (~ 30%) 

No y-rays could be identified from the previously reported 5 
level at 1426 keV (Bl 67). The production cross section of y-rays 
from this level to the 4 state of the g.s. band is, as we have 
seen in the case of Sm, small. The fact that the weight of the Sm 
sample was twice that of the Gd possibly explains why we were able 
to see the 5 » 4° transition in the case of Sm. 

The 1389 and 1215 keV y-rays were assigned to the decay of a 
level at 1464 keV to the 2* and 4* members of the g.s. band. This 
level has previously been observed in (p,p') and (d,d') studies 
and was assigned by the latter a value of J" = 3°. The observed 
decay is consistent with this assignment. 

The existence of levels at 1569, 1587 and 1591 keV was based 
on the observation of de-excitation y-rays from these levels to 


the 0° and 2° states of the g.s. band. 
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No quantitative analysis has been made for levels above 1400 
keV excitation due to the weakness of the observed y-rays. 
4.4.6 Level scheme and band structure of '*'Gd 

The decay scheme of *°®°Gd revealed from the present study is 
Shown in fig. 4.22, together with previous (p,p') and (d,d') 
results. The doubted levels at 946 and 1016 keV reported in the 
(d,d') work of Bloch et al. (Bl 67) were not seen in the present 
experiment. Except for the 5 level at 1426 keV observed in the 
(d,d') study, we have identified de-excitation y-rays from all the 
previously reported levels in the 900 to 1600 keV excitation region. 
In addition, seven new states have been proposed and included in 
the level scheme of *°°Gd. 

Analysis of the present results has confirmed the previous 
spin assignment of 2° and 3° for the levels at 989 and 1290 keV, 
respectively. Also, the spins of the 1058 and 1225 keV states 
have been determined to be eH and 1 , respectively. The new 
level at 1351 keV has been assigned a J" = 1°. In addition, 
tentative spin assignments have been made for the other levels shown 
Mi tigs 4222. 

Rotational bands built on the y- and octupole vibrational 
modes were suggested by Bloch et al. (Bl 67). The present study 
has completed the information about members of these bands as 
indicated in fig. 4.22. No evidence of the B-band was found. 

This situation is similar to that reported in the case of !°°Gd 


(Sh 71). The 3. level at 1464 keV might be a K = 1 octupole 
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Figure 4.22 


Level decay scheme of '°°Gd. Members of 
octupole and y-bands are also shown. The 
previous (p,p') and (d,d') results are from 
Shelton (Sh 66) and Bloch et al. (Bl 67), 


respectively. 
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vibration. The 1 member of the associated rotational band could 


be either the 1351 or 1378 keV level. 


4.5 Interpretation of Results in Terms of the Collective Model 


4.5.1 Background 


According the the collective model of Bohr and Mottelson 
(Bo 53), nuclei whose equilibrium shape deviates from spherical 
symmetry can undergo two essentially different modes of excitation -- 
rotational and intrinsic. The former is associated with a collec- 
tive motion which affects the orientation in space; the latter may 
be associated with the excitation of individual particles or with 
collective vibrations of the nuclear surface. Two assumptions 
were made: 

(i) The nuclear shape is axially symmetric. Thus the 
rotational motion can be characterized by the quantum 
numbers I, K and M representing the angular momentum, 
its projection on the nuclear symmetry axis, and its 
projection on the space-fixed axis, respectively. 

(ii) The rotational and intrinsic motions of the nucleus do 
not disturb each other. This adiabatic assumption 
corresponds to the existence of an approximate nuclear 
wave function given by a simple product of the intrinsic 
and rotational wave functions. 

The states in a rotational band are characterized by the same 


intrinsic wave function and are labelled by different values of I. 
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In an even-even nucleus, the ground state has K = 0 and symmetriza- 


tion of the wavefunction limits the rotational band to 


igs Telco [4.2] 


where Jf is the moment of inertia that depends on the nuclear 
deformation. 

Vibrations of axially symmetric deformed nuclei may be 
characterized by the quantum number A corresponding to the 
multipole order for small nuclear eccentricities [the parity 
of the vibration is Casual and by v, corresponding to the 
component of vibrational angular momentum around the symmetry 
axis. The lowest lying vibrations of a deformed axially 
symmetric nucleus are quadrupole vibrations, \ = 2 (g-vibrations) 
and one with y = 2 (y-vibrations). \y = + 1 can be shown to be 
the same as a rotation. 

Thus the g-vibration, having K = 0, is an axially symmetric 
vibration in which the eccentricity of the ellipsoid of revolution 
changes. In the y-vibration, the axial symmetry is destroyed, 

K = 2; that is if one looks along the symmetry axis, the outline 


is an ellipse. 
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In deformed nuclei the lowest lying negative parity vibrations 
should be octupole vibrations (\ = 3, v= 0, +1, + 2, + 3). On 
the other hand, for a spherical nucleus the strength of this 
vibrational mode should be concentrated in a single I" = 3° 
state. Octupole vibrations with K=0andI=1,3,5, etc. 
have been observed in many deformed nuclei. 

Deformed spheroidal nuclei can simultaneously rotate and 
vibrate. Thus, for each vibrational mode there is a band of 
rotational states whose energies are given by 


Bid a Peery ee | [4.3] 


since K is the z-component of I, I must be equal to or greater 
than K. The lowest energy state is then the one with I = K. 
For y-vibrations K ="2, so the values of [° =2,3 54, ... 
and for g-vibrations K = 0, where the values of I” are 0, “hs 
4”, (shite 

The adiabatic assumption of Bohr and Mottelson allows the 
reduced transition probability of a given multipole radiation L 
between two states to be expressed as a geometrical factor 
depending only on the angular momenta I, K and L, and a matrix 
element involving integrations over the intrinsic wave function 
of the initial and final states. When one compares the reduced 


transition probability from a state i to different members f, f', 


. of a rotational family, the factor involving the intrinsic 
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wave functions is the same. One thus obtains (Al 55) 


2 
B(L, I,K. * I Ke) ; <I.LK; Ke-K, [I Ke> 


so bh OARS Lia miele alone 
B(L, I,K; > IpsKe) <ILKs Ke-K, |Tpike> 


where < >is the Clebsch- Gordan coefficient for the addition 
of angular momenta I. and L to form the resultant I. This 

relation, of course, also holds where the states i and f, f' 
belong to the same rotational band. 

The simple separation of intrinsic excitations and 
collective rotations can be realized in the limit of large 
deformations, where the rotational motion is so slow that it 
does not disturb the nucleonic configuration or distort the 
nuclear shape. With decreasing deformation and increasing 


rotational frequency, the intrinsic structure is excited by 


the rotation motion, and the intrinsic quantum numbers are no 


longer constants of motion. This implies a modification in 


both the rotational spectrum (4.1) and the transition probab- 
ility ratios (4.3). The correction term for rotational ener- 


gies is usually taken to be proportional to I* (1 + 1)?, as is 


characteristic of the rotation-vibration interaction in 


molecules. The modified B(E2) for a transition from, say, the 


B-band to the ground state band can be written as, to the first 


order in spin-dependence (Fr 69) 
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[4.4] 
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B(E2sK=0,1, > K=0,1,) = <1020/ 1,0 |ma®}2 (1 4z, 


[Ig(1,#1) = Ip (Ig+1)] 3? [4.5] 


where MP 


is the leading-order matrix element that connects the 
B-band with the ground state band, and Zg is called the band 
mixing parameter. A similar expression can be derived for the 
y-band mixing with the ground state. Assuming that the intrinsic 
quadrupole moments for bands are equal, then Z. and 4y can be 
determined from a ratio between relative B(E2) values of two 
transitions depopulating a particular B- or y-band level. The 
magnitude of these mixing parameters provides a measure of the 
adequacy of the adiabatic assumption. 

4.5.2 Comparison of experimental and theoretical B(EL) ratios 

Using the branching ratios listed in tables 4.2 and 4.3, a 
series of transition probability ratios were calculated for y-rays 
connecting vibrational bands with the ground state rotational band. 
For transitions which might contain an ML component, it was assumed 
that it is small and can be neglected. This is reasonable since 
the transition probability of the competing EL is large. 

The experimental B(E1) ratios for y-rays depopulating the 
octupole band observed in *°*Sm and *°®°Gd are compared in table 4.4 
with the predictions of the adiabatic intensity rule (eq. 4.4). 
Good agreement is obtained for a K = 0 assignment. 

Table 4.5 shows the results for the y-band. The mixing para- 
meter Z,, was obtained using eq. (4.5). Bes et al. (Be 65) have 


calculated Z, for 15*Sm and '©°Gd to be 0.04 and 0.025, respectively. 
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TABLE 4.4 


Experimental and theoretical ratios of transition probabilities 
from members of the octupole band in '°*Sm and '°°Gd 


an BiGhalg > 14)/B(EVs1, > 14) 
O 
Nucleus : 
on et Experiment 
2 A er RON Theory* 
This work D'Auria et al. 
oe 
Loe Sin ——_*, 1.83 + 0.1 1.96 1.997 
3.374 
Ses 
2 1.19 + 0.06 1.25 ie33 
ce 
a ot 
160Gq es 2.07 + 0.1 1.997 
S + ee 
3° > 4 
peg 
2 1.32 + 0.07 132 
3° +2 


a) Predictions from adiabatic symmetric-rotor model (eq. 4.4). 
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TABLE 4.5 


Experimental and theoretical ratios of transition probabilities 


from members of the y-band in !°*Sm and 1°°Gd 


ae BUE Oat) /E(b2s Leena 
bial lI al ce ee one Abo ee ot 2 
Nucleus Te ce ANG) 
|e rece te ¥ 
Y g Experiment 
Theory? 
This work D'Auria et al. 
rie 
+S © 3, 0.50 + 0.04 0.4 0.4 4.0 + 0.32 
sew 
ees 
182 Oell Hs06 1.43 4 OU Ono 
oO 
> qt 
160Gq 0.54 + 0.03 0.4 5.4 + 0.33 
> at joy. 
> oF 
Z 1.75 # O01 aw Fuad 0220 
"alee 8| 
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Ei Ae Agee Ons erly 
Le ~ 


a) Predictions from adiabatic symmetric-rotor model (eq. 4.4). 


a ; on on 5 ae - vt , Maran 
yen a Ri Pa ‘ r y a nie Fae | : a) 
iy i} aye yy! f ey i - mS t 
Wa i , wy” , 
, a bi Ae 7 
— { - F 
F . } . - i - : :. j crs &. fs TBAT 


aeisti Fdsdova. noid tense 40 otter (9 
ae aan ‘hers ie aie rt ey ae to 2 


ae 


SU) 


Tey ER 


iris } E eit ; . xis . * op aaa a mit ; 7 . 
oe See ee (et i a J er Mn | 
*0F mM va i | nn F mace ete . - “ Le See | Jule a % m4 ORES 2 ~ a — 
bine | -Snanbyaqy a oe ie 


Synositt 


: ts 25 Peels. se ean al | 


C.04 0.8 “ £6 fo 


0+ 0.8 fA a oa | 


: ; + ‘ ~ oe 
Ce a a ane fee ; 


O80 4 be ear 


. od 4 | 
; f A ; ie 5he oT j 
eee es 


. 


i f8-® pd) rabon abvoa<3? ag 


9] 


TABLE 4.6 
Experimental and theoretical ratios of reduced E2 transition 


probabilities from members of the B-band in /°*Sm 


Bie 2 ale aba taal) 


ea 1s. 58 gooe, Dereeeens wren the ac sabett 
oe Tak Experiment Z, x 10 
B g a Theory® 
This work D'Auria et al. 
oa 2.134 0.16 2.07 1.8 0.63 + 0.05 
Cee = 
nee 0.45 + 0.03 0.47 0.7 Sue 
2>2 7 
aa 4.74 + 0.33 4.2 2.6 aes | 


AVES) 180 


a) Predictions from adiabatic symmetric-rotor model (eq. 4.4). 
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The present values are in good agreement with these predictions. 


The B(E2) ratios from members of the B-band in '°*Sm are given 
in table 4.6. There is an apparent inconsistency with the angular 
momentum dependency of B(E2) ratios, which is characteristic 
of B-bands of nuclei in this region (Ri 69, Fr 69). 


68 
for ***Sm (Ri 69, Ru 71). This is expected since the latter 


The values of Z, and z., obtained for '°*Sm are < half those 


nucleus is less deformed. Fraser et al. (Fr 69) and very recently 
Diamond et al. (Di 71) have determined the B(E2) values for 
transitions within the ground state band of '°*Sm. They 
found that it is in good agreement with the adiabatic values and 
if there is mixing it is very small, thus, supporting the present 
results. 

The smallness of band mixing parameters in ‘°"*Sm and 
160Gq indicates that these nuclei can be represented by 


rigid rotors. 
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CHAPTER V 


CONCLUSIONS 


The small sample method that has been developed in the present 
study has proven to be very useful for the investigation of the 
(n,n'y) reaction when only small quantities of the scattering 
material are available. It has been shown that the physical 
arrangement of source and scatterer used in this method performs 
a good energy average for inelastic cross sections, which can, 
therefore, be directly compared to the predictions of the statistical 
model. 

Application of the method to the study of the structure of 
15%*Sm and 1©°Gd has yielded information about excitation energies, 
principal decay modes and inelastic cross sections of levels in 
these nuclei. Interpretation of excitation functions in terms of 
the statistical model, together with the observed decay modes, has 
confirmed or determined level spins. The results are summarized 
in figs. 4.11 and 4.22. 

A spherical optical model potential was found to be adequate 
in predicting inelastic cross sections on these highly deformed 
nuclei. 

Previous reaction studies preferentially excited different 
levels in these nuclei. The fact that most of these levels were 


observed in the present experiment reflects the usefulness of the 
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‘CN' (n,n') reaction in populating levels independent of their 
structure (subject, however, to energy, angular momentum and 
parity conservation restrictions). 

Transition probability ratios obtained from the present study 
indicate that any coupling between the vibrational and rotational 
modes of excitation is very small. This means that the !°"Sm 
and '©°Gd nuclei are good rotors. 

An improvement of the experimental method described in this 
thesis could be made. As mentioned in Chapter IV, the low-energy 
region of the y-ray spectrum (below 500 keV) was badly masked by 
the intense y-rays associated with the ’Li + p reaction that was 
used aS a neutron source. By switching to the T(p,n) reaction 
it would be possible to obtain spectra free from these background 
y-rays. Identification of low-energy y-rays would be of assistance 
in assigning interaband transitions and accounting for cascade 
feeding. However, other background problems would arise from the 
use of a tritium gas cell. For example, one would probably see 


neutron and proton indiced y-rays from the gas cell material. 
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MEASUREMENTS OF NEUTRON ANGULAR DISTRIBUTIONS 
FROM THE 7Li(p,n)7Be REACTION™ 


S.A. Elbakr, I.J. van Heerden*, W.J. McDonald 


and G.C. Netlson 


Nuclear Research Centre 
The University of Alberta 
Edmonton, Alberta, Canada 


The absolute zero-degree differential cross sections of the 
7Li(p,no) 7Be and the “Li (p,ni) “Be reactions have been measured 
with a time-of-flight spectrometer at 0.1 MeV intervals from 
2.20 to 5.50 MeV. Angular distributions of neutrons from both 
reactions were also measured at 0.2 MeV intervals from 2.60 to 
5.40 MeV for the no neutron group and from 3.20 to 5.40 MeV for 
the ni group. Absolute differential cross sections were obtained 
by normalization to the zero-degree measurement. The centre-of- 
mass distributions were fitted with a Legendre polynomial series. 
Total cross sections were obtained by integrating the angular 
distributions. The comparison between present and previous 
results has revealed some disagreement in the shapes of angular 
distributions and in the magnitude of total cross sections for 


the ’Li(p,ni)’Be reaction. 
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1. Introduction 


The 7Li(p,n)7Be reaction has been used in the University of Alberta 
laboratory to study (n,n'y) reactions using small scattering samples in 
a close geometry!). Analysis of the y-ray data requires an accurate 
knowledge of the neutron angular distributions. The published informa- 
tion on the neutron angular distributions2-3) was not sufficiently 
complete or accurate enough for the analysis of our (n,n'y) data. To 
generate a complete set of angular distributions, a systematic 
measurement of the absolute differential cross section as a function 
of both energy and angle for the ’Li(p,n,)7’Be and 7Li(p,n,)7Be reactions 
has been carried out at proton energies from 2.60 to 5.40 MeV and 
laboratory angles between 0° and 150°. The measured angular distribu- 
tions were fitted with a Legendre polynomial series yielding a set of 
coefficients at the different proton energies to be used in calculating 


differential cross sections at any angle in this energy range. 


2. Experimental procedure 


The measurements were made with neutron time-of-flight spectro- 
meters using pulsed proton beams from the 6 MV Van de Graaff acceler- 
ators of the University of Alberta and of the Southern Universities 
Nuclear Institute (SUNI). 

At Alberta, the absolute differential cross sections of the 
7Li(p,no)7Be and 7Li(p,n,)7Be* reactions were measured at zero degree 


for proton energies from 2.20 to 5.50 MeV at 0.1 MeV intervals. 
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Angular distribution measurements at 20° intervals were then made at 
proton energies from 2.60 to 5.40 MeV for the ng neutron group and 
from 3.20 to 5.40 MeV for the n, neutron group. 

The neutron time-of-flight system used has been previously 
described). Top terminal pulsing and a Mobley bunching magnet 
provided proton pulses of < 0.5 ns. The neutron detector consisted 
of a cylindrical glass cell 9.5 cm diameter x 1.9 cm thick filled 
with NE 213 liquid scintillator viewed by an RCA 4522 photo-multiplier, 
and was placed at a distance of 3.473 m from the target. General 
background was reduced by pulse shape discrimination and the overall 
time resolution of the spectrometer was about 1 ns. 

The neutron detection efficiency was calculated using a Monte Carlo 
program written by N.R. Stanton®). This program was adapted for liquid 
scintillators by incorporating the data of Smith et al.°) for the light 
output from proton recoils in the NE 213 liquid scintillator. For each 
neutron energy the program gives the efficiency as a function of 
neutron detection threshold which was independently determined. The 
results of Chastel et al.’) have shown that the calculated efficiencies 
are in good agreement with measured values. 

The lithium target was made by evaporating LiF on a 0.01 cm 
tantalum backing to a thickness of approximately 200 yg/cm*. The 
total thickness has been checked by measuring the width (FWHM) of a 
low energy neutron peak and subtracting in quadrature the time 


resolution of the spectrometer. The tables of Janni®) for proton 
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energy loss were then used to determine the target thickness in mg/cm2. 
This gave a value agreeing to within 10% with that determined by 
weighing. 

To ensure that the proton beam was positioned on the same target 
spot throughout the measurements, the beam was passed through a pre- 
target slit system. This consisted of two gold plates independently 
mounted on two micrometers and connected to two microammeters. A 
negative bias of 300V applied to each plate prevented secondary 
electrons from leaving the target. The Mobley magnet current was 
adjusted during the experiment to ensure a balanced minimum current 
On the plates. 

For the excitation curve measurements a beam current integrator 
was used as a controlling monitor. The target was air-cooled. To 
check the internal consistency of the results and whether any target 
deterioration had taken place, various cross section determinations 
were repeated at randomly selected proton energies. No significant 
differences were observed. 

For the angular distribution measurements a second time-of-flight 
system was used as a monitor. The neutron detector for the monitor 
system was positioned at a fixed angle of -30° to the incident proton 
beam. It was found that for each angular distribution measurement, 
the ratio of monitor counts to total electric charge on target 
remained constant to within the statistics. This indicated that the 


target was probably uniform, and that the thickness determination by 


lit 


¥} 


a 
- 
— 
7S 
7 
= 
et 


1 =i} a 7 
) uh 
, 7 ' 
; = 
~ ; 7 
- 6 
» ‘ ; 
Be: 
iy 


. Sa \om i deena pit ‘Jopted bang en tnrrad ab a beau agi 


ud brett widteb tony hws el ntrd tw vt ons i 6 pen 


i Vi rf ¢ 
2 bie * . 
: : t The : 4 ie bios ae allie 
_ Septet sms2-snd no “bsaotsteaon 2aw: med Midi ti Sail. oe Ni 
ee as rain | * ' ia 


. : a 
: ‘ 7 } , "3 - ef 4 bs er o 
“Bq 6 TOMO “bee2nq ew mpad ons. ediems ywesam ors suodpudtdd oge 
: a 
visnebne Shabir f 28 f6 iq blop owe 79 Desei Lene ‘Bat Kanere eye ste! 7“ 


~ A  .2astemispoyotm ows. oF bstoanno ins sat iio owt 0 bedi 


Viebnosse besnove 1g stefq does oF bo} fags: Wooe to: este vs 


~ 
ay. 
ay 


GRSTND ‘Sonpasm YSTdoM SAT sprss” ‘of itwost mort sd 


toetyuoMunMinim bsonelsd a sywehe oF anise bAt ‘enfrub ane r 
| : jiae< 


sal 

‘ eo } a. 
‘ . mh” | zeae he, 9 wire 
e il a a 
Pe A, s ae er 


rods ipso f tayo msod 5 2tnom vesom sv13) otbRd faxo ont 40} 


Vise 


‘ingle & 26 dere 


OT bal 003+i6° 26W 39 ete) 1 | — ott 


= lte 


2iotteninatab word 22049. 2UOMIBY . 898 re 6 
a aT? i sade iat 
sue zag ai notord bstosfee. wie 


oy 


eo iat a ; 
norlt-to-emtt b 10992, 6 einelEsGw not stun ot 10° fit - 
vot rnon sdt yor vodasteb! ponsusa odT 03 and 28h 


one.) 


sp a 


flodova tnebtont ori os “Oe- to- ‘fens bax 6. 


ethomeiue2samn notdudnazb, eer at 
; ie Wray E i 
-t$py6s no avari sivas fatod-o ‘= Ve al t 

4am | a | 

edt tent bassoibnt ahaT’ ‘aotta iad Site 


wv 


yd er ee eeomiot ht ont dat 
} ane | 


Ja: 


weighing was justifiable. Relative cross sections were determined by 
normalizing the yield in the main detector, corrected for efficiency, 
to that in the monitor. 

At SUNI the angular distributions were remeasured using a similar 

neutron time-of-flight spectrometer at the same proton energies as 
before. Top terminal pulsing and a klystron bunching system provided 
proton pulses of ~1.0 ns. The main neutron detector was a cylindrical 
glass cell 12 cm diameter x 4 cm thick filled with NE 213 liquid 
scintillator mounted on a 58 AVP photo-multiplier and it was placed at 
a distance of 4.08 m from the target. Pulse shape discrimination 
reduced the general y-ray background. A standard Hansen-McKibben long 
counter with known efficiency?) placed at 90° was used as monitor. 
The relative efficiency of the neutron detector was measured at low 
energies (threshold to 3.5 MeV) using neutrons from the 7Li(p,n)’Be 
reaction and direct comparison!9-11) with the known response of the 
Standard long counter. 

It was found that the two sets of angular distribution measurements 
agree within statistics. They were then normalized to the zero-degree 


absolute differential cross sections. 


3. Results and Discussion 
Fig. 1 presents the zero-degree differential cross sections for 
both the ground state and first excited state groups. Statistical 


errors are < 3%, while the error on absolute values is + 12%. The 
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latter arises mainly from uncertainties in target thickness and detector 
efficiency determinations. The present results agree, within error, 
with the absolute values of Gabbard et al.12) for the zero-degree 
differential cross sections of all neutrons. They are also in fair 
agreement with the results of the relative cross section measurements 
of Bevington et al.2) and Borchers and Poppe?). 

The centre-of-mass angular distributions were fitted, using a 


least squares program, to a legendre polynomial series of the form: 


o (6) = } A. Pi (cosé) ey 


Terms up to and including P, were found sufficient to yield a good 
fit. Figs. 2 and 3 show the results where the smooth curves through 
the data are the least squares fits. Statistical errors on the data 
points are < 5%. The coefficients resulting from the fits are given 
in Table 1. 

The (p,n,) angular distributions vary in a smooth manner from 
forward peaked at low energies to symmetric around 90° at energies 
exceeding 5.0 MeV. The (p,n,) distributions tend to be peaked in a 
backward direction. In order to compare the shapes of the present 
angular distributions with those of Bevington et al.2) and Borchers 
and Poppe3), we have fitted their results with our least squares 
program. The ratios A /Ag for both neutron groups are plotted as 


a function of proton energy in figs. 4 and 5. It is clear that 
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there is some disagreement between the present and previous results 
regarding the shapes of the angular distributions. 
By integrating the Legendre series [1] over the solid angle, the 


total cross section is given by 


Clee Aaa [2] 


Total cross sections obtained using this equation are shown in 

fig. 6, where they are also compared with other measurements. The 
present measurements of the ground state total cross sections agree 
with the previous measurements of Bevington et al. and Borchers and 
Poppe, but for the first excited state the present total cross 
sections are about 20 -30% smaller. Presser and Bass!3) have also 
recently determined the total cross section for the (p,n,) reaction 
by measuring the 90° yield of the 0.431 MeV y-ray. The present 
results are 13-22% smaller than their values. 

In fig. 6 we have also compared the total cross sections for 
reactions leading to the ground and first excited states of ’Be with 
the absolute measurements of Gibbons and Macklin!*). Their data were 
obtained by thermalization of neutrons in a large graphite sphere for 
which the absolute neutron efficiency was very accurately measured. 
The absolute cross sections were based on weighing a thick LiF target, 
380 wg/cm2, and were accurate to + 5%. The present total cross 


sections are accurate to + 12%, and in general are smaller than those 
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of Gibbons and Macklin'*) by < 10% below 5 MeV, well within the 
uncertainties of the absolute cross section scales. 

It is obvious that the ratios I,/Ig of the intensities for 
the (p,n,) and (p,ng) total cross sections obtained in the present 
work are smaller than those obtained by Bevington et al.*) and by 
Borchers and Poppe?). The ratios obtained at 0° in the present 
work and those obtained by the above authors and also by Cranberg!*) 
are in good agreement with one another. Differences in the angular 
distributions obtained must therefore account for the differences 
in the I,/I,g ratios for the total cross sections. The present 
data is based on two independent measurements of the same set of 
angular distributions. Although different methods were used to 
determine the relative efficiencies of the neutron detectors, the 
two sets of data were consistent within the errors. Furthermore, 
the present measurements were carried out with high resolution, low 
background time-of-flight spectrometers. We therefore feel that 
our results are more reliable than those of Bevington et al.*) and 


Borchers and Poppe?). 
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Figure Captions 


Figure 1 


Figure 2 


Figure 3 


Figure 4 


Figure 5 


Figure 6 


Absolute zero-degree differential cross sections. The 
upper and lower curves are for neutrons leading to the 
ground and first excited states of ’Be, respectively. 
Absolute errors are represented by bars for the no 


neutron group and by the size of points for the n, group. 


Centre-of-mass differential cross sections for the 
7Li(p,no)’Be reaction. Each distribution is offset by a 
certain amount, so that the amount indicated at the right 
must be subtracted from the ordinate to obtain the 


correct magnitude. 


Centre-of-mass differential cross sections for the 


7Li(p,ni)’Be reaction. (See also caption of fig. 2). 


Legendre coefficients for the centre-of-mass angular 


distributions of the ’Li(p,no)’Be reaction. 


Legendre coefficients for the centre-of-mass angular 


distribution of the 7Li(p,n,)’Be reaction. 


Total reaction cross sections. The data of Gibbons and 
Macklin?!) were obtained by thermalization of neutrons 
in a large graphite sphere. All other data were obtained 


by integrating angular distribution of neutrons. 
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Geometrical Integration of Neutron Flux 


In chapter III, the average production cross section has 


been defined as (eq. 3.6) 
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Apart from the exponential term, the integration of (B.1) can 


be written as 
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The limits of integration are those given in chapter III. 
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when z = 0 cot a 8 = a 
and when z = T = op cot 6, 6, = tan” o/T 
Eq. (B.3) can then be written as 
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Thus the integration of (B.1) is given by 
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The integration (B.4) can be written explicitly as: 
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when the limits of integration are applied, we get f(p,T) to be 
used in eq. (B.5). The integrations over p,d must be carried 


out numerically. 
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